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We study the steady state of a zero-temperature Bose gas near a Feshbach or photoassociation res- 
onance using a two-channel mean-field model that incorporates atomic and molecular condensates, 
as well as correlated atom pairs originating from dissociation of molecules into pairs of atoms. We 
start from a many-body Hamiltonian for atom-molecule conversion, and derive the time dependent 
version of the mean-field theory. The stationary solution of the time dependent model is rendered 
unique with an approximation that entails that all noncondensate atoms are correlated, as if emerg- 
ing from dissociation of molecules. The steady state is solved numerically, but limiting cases are 
also found analytically. The system has a phase transition in which the atomic condensate emerges 
in a nonanalytic fashion. We quantify the scaling of the observable quantities, such as fractions of 
atomic and molecular condensates, with the detuning and the atom-molecule conversion strength. 
Qualitatively, the dependence on detuning rounds out with increasing coupling strength. A study 
of the thermodynamics shows that the pressure of the atom-molecule system is negative, even on 
the molecule side of the resonance. This indicates the possibility of mechanical instability. 

PACS numbers: 



I. INTRODUCTION 

Photoassociation [TJ [5] and especially the mathemat- 
ically equivalent Feshbach resonance [3] have been 
main themes in the physics of quantum degenerate Bose 
and Fermi gases for a while. The earliest experimental 
work was carried out on bosons. The widest-known ex- 
amples are probably Feshbach resonance in a BEC [5J, 
two- photon photoassociation [6] in a BEC, purported [7^ 
mechanical collapse of a BEC when the atom-atom in- 
teraction is made attractive by using a Feshbach reso- 
nance [5], and Ramsey fringes in the conversion between 
atomic and molecular condensates |9]. 

In an aspect that is often concealed by phrases to 
the effect that the Feshbach resonance is used to tune 
atom-atom interactions, at least in principle photoasso- 
ciation and Feshbach resonance always involve conversion 
of atom pairs to the corresponding diatomic molecules. 
These molecules are in a highly excited vibrational state 
and are prone to coUisional quenching. In a Bose gas 
the molecules tend to be short-lived (~ 10 ms). How- 
ever, it turns out [101 E] that the diatomic molecules 
created at the 834 G Feshbach resonance from fermionic 
^Li atoms may persist for seconds. Thermal-equilibrium 
experiments are feasible, which has spurred an enormous 
experimental and theoretical interest. Molecular conden- 
sates are now prepared routinely [121 US] • Of particular 
interest to our theme of atom-molecule conversion is a 
measurement of the equilibrium fraction of molecules as 
a function of the magnetic field [14| . while a compari- 
son of the thermodynamical properties between exper- 
iments |15j and quantum Monte Carlo simulations |16j 
serves as an illustration of the present state of the art in 
the modeling of fermion systems. 

Given that the initial Feshbach resonance experiments 
were time dependent in an essential way, the correspond- 



ing theories were time dependent as well. Our approach is 
based on field-theoretical modeling of atom-molecule sys- 
tems [T71 HSl US] , but there are methods springing, e.g., 
from theory of molecular structure, as reviewed recently 
in Ref. [20]. For instance, a good agreement with the 
atom-molecule Ramsey fringe experiments in a BEC [S] 
has been reported by several groups jUJ |^ |^ [511 [IS] . 
Our time dependent method has produced [26] a passable 
theoretical description of atom-molecule conversion in an 
experiment (27j , and yielded valid |28j predictions about 
the temperature dependence and the maximum value of 
the conversion efficiency also for fermionic atoms. 

With the emergence of thermal-equilibrium experi- 
ments in Fermi systems, the equilibrium-oriented theo- 
retical machinery of condensed matter physics has been 
brought to bear [T6j [29l [30] . On the other hand, we have 
noticed [31 that the stationary solution of our time de- 
pendent formalism for fermions [26 may also serve as a 
zero-temperature thermal equilibrium theory. The sta- 
tionary solution is not unique, but it may be made so 
by assuming that all fermions appear in correlated pairs, 
as if from dissociation of molecules. The result turns 
out to be a variant of the atom-molecule version of the 
BCS theory [32l [33||34l |35] for Fermi gases. We [31^ and 
others [501 [SS] have reported favorable theoretical com- 
parisons with the equilibrium fraction of molecules in the 
experiments |14j in which the magnetic field was varied in 
the neighborhood of the broad ^Li Feshbach resonance. 

The immediate purpose of the present paper is to in- 
troduce a similar "BCS" theory for bosonic atoms at 
zero-temperature thermal equilibrium. Basically, we take 
our time dependent theory ^22i [25l [37] and find the 
time-independent solution. As with fermions, the sta- 
tionary solution may be rendered unique with a pairing 
approximation, which, however, turns out to be more 
subtle for bosons than for fermions. As the secondary 
goal, we discuss the technicalities of our work on both 
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bosons |221 [251 EZ] (^''^d fermions |3T] that were not 
detailed in the original letter format publications. 

The main qualitative finding is that the atom-molecule 
system exhibits a phase transition in which the atomic 
BEC emerges in a nonanalytic fashion when the detun- 
ing, the atom-molecule energy difference, is varied. In 
the limit of weak atom-molecule conversion, for instance 
in the limit of a very dilute gas, the phase transition 
is at the position of the two-body Feshbach resonance, 
but for increasing atom-molecule conversion strength it 
moves to the molecule side of the Feshbach resonance or 
photoassociation. We characterize the phase transition 
and the variation of quantities such as the molecule frac- 
tion with both the detuning and the atom-molecule con- 
version strength. Broadly speaking, the dependence on 
detuning rounds out with increasing interaction strength. 
Finally, we investigate the thermodynamics of the atom- 
molecule system. Our present model with a Feshbach or 
photoassociation resonance but no attendant background 
scattering length produces a puzzling surprise: The cal- 
culated pressure of the gas is negative for all parameter 
values, which indicates mechanical instability. 

II. FORMULATION OF THE PROBLEM 
A. Hamiltonian 

We model conversion of bosonic atoms into bosonic 
diatomic molecules using the momentum (wavevector) 
representation Hamiltonian [111 ISH] 

y = [ck ajak + {6+ |£k)&k^k 

k 
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as appropriate for a free (non-trapped) gas. Here at 
and 6k are boson annihilation operators for atoms and 
molecules with momentum hk., and fiet = h^]i.^/2m is 
the kinetic energy for an atom with wave vector k. For 
a molecule with twice the mass of an atom, the kinetic 
energy for a given momentum is half of the energy of 
an atom. The detuning S gives the energy difference be- 
tween a stationary molecule (k = 0) and two stationary 
atoms in the form M. This parameter is varied in a 
Feshbach resonance by varying the magnetic field, and in 
photoassociation by tuning the frequency (frequencies) of 
the laser(s). 

Atom-molecule conversion entails that a pair of atoms 
is either converted to a molecule or a molecule is dissoci- 
ated to a pair of atom, all the while conserving the mo- 
mentum. The governing matrix elements are denoted by 
Kp^q. Ordinarily we deal with s-wave processes that dom- 
inate at low temperature. By virtue of translational and 
rotational invariance, the corresponding coupling matrix 
element is then of the form Kp q = k (|p - q|). Fur- 
thermore, we mostly write the coupling matrix element 



as a once-and-for all constant Kp.q = k. This means 
that atom-molecule conversion is a zero-range contact 
interaction; in the position representation for the atomic 
and molecular fields </> and ijj, the corresponding term in 
the Hamiltonian density would be oc K['i/'^(r)^(r)0(r) -|- 
(f)^r)(f)^r)^p{r)]. The contact interaction model is conve- 
nient for our aims for two reasons: it reduces the number 
of parameters to consider, and endowing the coupling co- 
efficients with a realistic momentum dependence would 
much complicate the solution of the model. The down- 
side is an ultraviolet divergence [37j , which will require a 
renormalization . 

Consider two atoms that may combine into a molecule 
in the center-of-mass frame. Then a molecular bound 
state is coupled to the continuum of the relative motion 
of the two atoms. By redefining the global phases of 
the continuum wave functions, one may always arrange 
things so that the coupling coefficients q and k are 
real, and k is also non-negative. This is assumed to be 
the case below. We write sums over the wave vectors in 
Eq. Q , which presumes box normalization of the dissoci- 
ated states. The quantization volume V therefore enters 
the matrix element, Kp.q oc 

Conspicuously missing from the model is a background 
scattering length that would prevail, say, far away from 
the Feshbach resonance. It could, and given the nature of 
our results, maybe should be included. But again, since 
our focus is on the nature of the system in the vicinity of 
the resonance (5 = 0, we keep the model and the attendant 
technical complications to the bare minimum. 

A break-up of a molecule produces two atoms, and 
it takes precisely two atoms to make a molecule. The 
Hamiltonian ([ij correspondingly has the conserved quan- 
tity 

^-E(«k«k + 25[6k), (2) 

k 

number of atoms plus twice the number of molecules. 
We use the value N of the invariant N to characterize 
the number of particles in the system, and occasionally, 
slightly inaccurately, call it the atom number. The total 
momentum 

P = hY,HaW + blbu) (3) 

k 

is likewise a constant of the motion. 



B. The dressed molecule 

Before proceeding to many-body systems, we investi- 
gate the solvable model with the invariant atom num- 
ber equal to two. Without restricting the generality, we 
also assume that the conserved center-of-mass momen- 
tum equals zero. The state space is then spanned by 
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vectors of the form [3T] 



E 



f3bl |0) 



(4) 



where v4(k) [with A{k) = k)] and (3 are complex 
numbers, and |0) is the particle vacuum. The time depen- 
dent Schrodinger equation follows from Hamiltonian ([T]) 
in the form 



iA{k) = 2ekA(k) - k(3 . 



(5) 
(6) 



If rotational invariance prevails at the initial time so 
that A(k) is only a function of |k|, or equivalently, a 
function of ek = e, the same symmetry holds at all later 
times. Moreover, in a two-atom system we may replace 
the sum over k by a continuum approximation without 
running into problems with the atomic BEC. We write 
the continuum approximation as 
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(7) 

We prefer quantities with the dimension of frequency over 
h times the same quantities with the dimension of energy, 
so that the integral runs over frequencies. Although it 
is silly if not misleading in the present case of only two 
atoms (N = 2), for future use we have defined the energy 
hep by 



h /67riV\ 
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It equals the Fermi energy for a single-component gas 
with the density N/V, but obviously has nothing to do 
with any physical Fermi energy. Instead, e^? is a measure 
of the density of the gas; the essentially unique frequency 
that can be construct out of density (N/V) for a quantum 
mechanical (h) gas of atoms (to). Below we will refer to 
ep without the h as Fermi energy. 

To cap our introduction of the notations, we define the 
frequency parameter characterizing the atom-molecule 
coupling 



n = Vnk. 



(9) 



It is essential to keep in mind that in this bound- 
continuum problem the analogy of the Rabi frequency 
f2 is proportional to the square root of density, f2 oc 
(N/V)^/'^. In line of what was said before of the co- 
efficient K, we also take > 0. The time dependent 
Schrodinger equation for the coefficients P and A{e(k)) 
finally reads 
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d€^/^A{e,t), (10) 
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The essence of the two-channel theory is to regard 
atoms and molecules as distinct though coupled degrees 
of freedom. The boson operators in Hamiltonian ^ cre- 
ate and annihilate bare atoms and molecules that would 
be observed if there were no atom-molecule coupling. As 
such, they would represent the observable atoms and 
molecules immediate after the atom-molecule coupling 
were switched off. In the case of photoassociation this 
could be achieved literally by switching off the lasers. 
For the Feshbach resonance an equivalent decomposition 
could be effected (in principle) by suddenly switching the 
magnetic field so far off the resonance that the atoms and 
the molecules effectively decouple. There are also exper- 
imental probes that directly see the bare molecules, for 
instance, by making use of optical transitions in the bare 
molecules [T3]. 

However, standard radio frequency spectroscopy at 
a Feshbach resonance |27l [39] probes transitions be- 
tween energy eigenstates of the system in the presence 
of the atom-molecule coupling, i.e., stationary states of 
Eqs. (10 1 and (11). These are superpositions of a bare 



molecule and pairs of bare atoms, and so we refer to the 
coupled system as the dressed molecule. 



RenormaUzation 



The energy eigenstates are obtained by inserting an 
Ansatz of the form (3{t) = e"*"*/?, A{e,t) = e~*'^*A(e) 



into Eqs. (10 1 and (111, which gives 



{lu - 2e)A{e) 
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n 



f3. 



deV~eA{e), (12) 



(13) 



Simple elimination of A{e) from Eq. (12 I using Eq. (13 1 
gives a relation to determine the eigenfrequency to, 
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de- 



Lj — irj — 2e 



(14) 



Here —rj, with r/ = 0+, is the usual imaginary part in the 
energy that needs to be added to handle the divergence 
of the integrand at 2e = uj. This practice is the same as 
if we took Fourier transformations of the time dependent 
equations and used them to study the evolution of the 
system forward in time. Such an asymmetry in the di- 
rection of time is not desirable if we are looking for true 
stationary states of the atom-molecule system. We out- 
line in Appendix [A| a method, following Ref. [10], to find 
the proper stationary states, but here the main issue is 



the ultraviolet divergence of the integral in (14) 



Physically, Eqs. (10) and (111 describe the coupling of 



the bare-molecule state to many (actually, a continuum 
of) atom-pair states. One obvious consequence is that, 
if there still is a bound state in the system, it is shifted 
in energy from the original bare molecular state. In the 
contact-interaction model the shift simply is infinite. 
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We renormalize as follows [3TJ 157], W e adopt an upper 
limit of the integral AI, write Eq. ( 14 1 as 



U! 



4e 
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M 



2e 
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oj ^ ill — 2e 
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(15) 



and let M ^ cx) at the end of the calculation. The right- 
hand side then converges nicely, but ostensibly not so 
the left-hand side; the detuning S gets modified by the 
infinite level shift. The idea of the renormalization is 
to incorporate the level shift into the definition of the 
energies, and take the renormalized detuning 
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(16) 



operators 



ia — —Kha\ 

ib = (56 — ^Koa — flkfl-k 



(18) 
(19) 

(20) 



Out of these primary equations one may form equations 
of motion for quadratic operator products, e.g.. 



^^(4ak) 



= K(6^aka_k - ' 



— 2ekaka_k 
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(21) 



(22) 



to have a finite value. 



In this way Eq. ( 14 1 turns into a well-behaved equation 



-S = 



. 3/2 , 



oo 2^ 

de ■ 



n oj + irj ~ 2€ 



(17) 



The salient point is that for (5 < Eq. (17 1 has precisely 
one real solution (with w < 0), and no real solution if 
^ > 0. The system has a true stationary state that does 
not evolve in time only for 5 < 0. In Appendix ^ we 
replace this statement with the more precise observation 
that the dressed molecule has a normalizahle stationary 
state if and only if J < 0. A bound state is found for 
the dressed molecule for 5 < 0, otherwise the dressed 
molecule only exists in a dissociated form as a pair of 
bare atoms with a component of the bare molecule mixed 
in. We take this to mean that the renormalized detuning 
^ = denotes the position of the Feshbach resonance in 
the two-atom system. 



We implement the mean-field approximation at this 
juncture by stating that in the equations of motion a 
and h are c- numbers not operators anymore [26' '3T . The 
quantum mechanical expectation values for the products 
then factorize as in this example, 



(aj^ak^) = (aj^ak)6 = (aj^ak)(6) 



(23) 



Applying the factorization to Eqs. (18l-(^22l gives a 



closed set of equations of motion involving the expec- 
tation values of the form (a), (6), (aj^ak), and (akfl-k)- 

For the convenience of the formulation we assume that 
the problem is effectively rotationally symmetric, so that, 
e.g., the expectation value (aka_k) only depends on the 
energy e = Ck- We write 



P(e) 



(OkOk) 



A{e) 



[aka-k) 



(24) 



C. Mean- field approximation 

We make the many-particle system solvable by a mean- 
field approximation, the idea of which is to treat possi- 
ble atomic and molecular condensates as classical fields 
not quantum fields anymore [ITl [HI [26]. Although it 
is not essential for the structure of the mean-field the- 
ory but rather a technical assumption to facilitate the 
analysis, we also assume that all molecules present in 
the system belong to a condensate of zero-momentum 
molecules. Hence, only the molecule operator bo = b 
is kept in the Hamiltonian. Moreover, to accommodate 
the corresponding atomic BEG, we already at this point 
track separately the zero-momentum atoms with oq = a. 
Using the Hamiltonian ([T]) we then find the Heisenberg 
picture equations of motion for the atomic and molecular 



It should be noted that P(e) stands for the expectation 
value of the number of atoms in a one-particle state with 
the energy he, not for a quantity such as the number 
of atoms per unit energy interval. We also define the 
amplitudes for atomic and molecular condensates 




so that jap and |/3p stand for the fractions of the atoms 
that are in the system as part of either the atomic or the 
molecular condensate. Finally, as the atomic BEG has al- 
ready been taken into account separately, the continuum 
approximation ^ should work as before. 

We finally have the equations of motion of our mean- 
field theory for atom-molecule conversion in a boson sys- 
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tern, 
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(26) 
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J deV~eA{e,t), (27) 



n 



2eA(e,t)--^[l + 2P(e,t)]/3(i), (28) 
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[P*{t)Aie,t)-P{t)A*ie,t)]. (29) 



The similarity to the notation we employed in the dis- 
cussion of the two-atom problem in Sec. |IIB| is no acci- 
dent. In fact, in the absence of a BEC of atoms, and 
assuming that the occupation numbers of the atomic 
states P{€,t) are negligible compared to unity, Eqs. (26) 
and (27) coincide with the two-atom theory Eqs. (10) 
and (111. We view the mean-field theory as the two- 



atom theory amended with the possibility of a BEC and 
Bose enhancement for the atoms. An analogous interpre- 
tation applies to the corresponding mean-field theory for 
fermions [31 . 

Writing the expectation value of the invariant atom 
number ^ in terms of the mean-field variables gives the 
equation 



2e 



3/2 



&yiP(e) = 1. 



(30) 



The left-hand side is indeed a constant of the motion 



under Eqs. (26l-(29l, so that the present mean-field ap- 



proximation successfully reflects a basic property of the 
Hamiltonian. Finally, one may write the expectation 
value of the Hamiltonian in the mean-field approxima- 
tion as 



hN 



2e 



3/2 



n 

2V2 
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2a . 3/3 



2e 



3/2 



de ^fe A{e) + c.c. . (31) 



Using Eqs. (26l-(^29^ it may be shown straightforwardly 



that, provided one is willing to subtract certain formally 
equal divergent integrals to obtain zero, the energy per 
particle he from Eq. (31 ) is also a constant of the motion. 



The divergent integrals are part of the issue of renormal- 
ization, which was already discussed above and will be 
revisited again shortly. 



We have made use of Eqs. (26l-(30) in different vari 



ations, notations, and approximations many a time in 
the past in our boson theories [T9j [22l [25l [37]. Simi- 
larly, we have repeatedly resorted to basically the same 
approach in the theory of the conversion between a 



two-species Fermi gas and the corresponding diatomic 
molecules [26l [31]. Besides the obvious absence of the 
fermionic condensate, the difference is that where the bo- 
son problem shows the factor l-|-2P(e) for Bose enhance- 
ment in Eq. ( |27[ ) , the corresponding fermion equation has 
the factor 1 — 2P(e) refiecting the exclusion principle. 



D. Pairing approximation for steady state 

Barring circumstances such as interference of the 
atomic BEC component of the atom-molecule system 
with another reference BEC, a multiplicative complex 
phase factor in the quantities a, (3 and A{e) is not ob- 
servable. Besides, it is obvious from the equations of 



motion (26l-(29) that a certain combination of exponen- 
tially evolving phases is self-sustained. Specifically, we 
search a stationary solution in the form |31j 



a{t) = e~'^*a, /3(t) 



f3, 



A{e,t)^Aie)e-'^^\Pie,t)^P{e), 



(32) 



where /Lt is a real frequency. It will turn out that is the 
chemical potential for the atoms in this system (and half 
of the chemical potential for the molecules) , but such an 
interpretation is not a given at this stage. In the rest 
of the paper we will again ignore the h, and call /i the 
chemical potential. 

Now, by a suitable choice of the zero of time we may 
always make the coefficient f] in Eqs. (32) real and non- 



negative, /? > 0; let us assume so from now on. To keep 
Eq. ( 29 ) valid with an Ansatz of the form ( 32 1 at all 
times is then only possible ii {(3 — or if) A(e) is real. 
Likewise, by Eq. (26), the amplitude a must be real. 



With these restrictions, the time independent coefficients 
must satisfy 



(3n 
71 
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ifi ~ e)Aie) 



- 0, 
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2e 
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J deV^A{e) 



[l + 2P(e)]/3, 



(33) 



(34) 
(35) 



and, of course, the norm condition (30 1. 



The unknowns in the steady state are fi, a, (3, A{e), 
and P(e). Thinking about a numerical solution with a 
discrete set of values for e, it is clear that there are many 
more unknowns than equations. The problem is the orig- 
inal Eq. ( 29 1 , which will not lead to any useful relation 
between A\e) and P(e) in the steady state. Additional 
conditions are needed to constrain the solution. 

The same dilemma came up in a system of two fermion 
species combining into bosonic molecules. We resolved 
it [31 by the assumption that the fermions only come in 
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pairs with opposite momenta and spins, as from dissocia- 
tion of molecules. This leads immediately to the relation 
between pairing amplitudes and occupation numbers 



1.01 



l(ckTC-ki)| = (Ck|Ckt) - (4^Ckt) , (36) 

or, in the notation of the mean- field theory in Ref. fST, 

\C{e)\'-[Pie)~P'{e)]^0. (37) 

C(e) is the fermion pairing amplitude analogous to A{e) 
of the present paper. With Eq. (37 1, the number of equa- 



tions was sufficient for a (presumably) unique solution. 
Moreover, while we have not mentioned this before, the 



left-hand side of Eq. (37 1 is a constant of the motion in 



our BCS style mean-field theory for fermions pTj . 

To address the corresponding boson case, let us take 
two momentum states ±, short for ±k, with the occu- 
pation numbers \n+n_). The most general completely 
paired state is of the form 



^Cn\n,n) , 



(38) 




FIG. 1: The ratio of the actual pairing amplitude A{x) and 
the quantity [P{x) + P'^{x)]^^^ derived from the occupation 
number of the boson state P{x) for the Poissonian paired 
state, (381 with (391, as a function of the real parameter of 

were exact, 



the state x = a. If the pairing approximation (|42 
this ratio would identically equal unity. 



with J2n = 1- Given the propensity of boson to 
Poissonian statistics, we take 



c„ = e 2 



(39) 



where a is a complex number. In the limit \a\ <^ 1 
this is a generic description for a situation when only the 
states |00) and |11) are occupied, and the latter with a 
much smaller probability. Similarly, in the limit of a real 
a = X ^ I we have a generic description of the state in 
which |c„| peak around n ~ x^, and c„ vary slowly as 
a function of n around the maximum. In fact, we cover 
the case |a| <C 1, too, if we just use a real and positive x 
in our argument, so that is how we proceed. Given the 
model, we have the expectation values 



P = (a^a+) = . 



(40) 

A EE (a+a_) = e x}_^ ^ . (41) 



Mindful of the sign differences between bosons and 
fermions, on the basis of Eq. (37 1 we expect a relation 
for bosons of the form 



\A{e)\'-[P{e)+P\e)] = 0. 



(42) 



To find out if it wor ks, we pl ot in Fig. [T] the ratio 
A/VP + P'^ from Eqs. (40 I and (41 1 as a function of the 
variable x. The maximum deviation of this ratio from 



unity is about three per cent. We surmise that (42 1 is 



a reasonable approximation between pairing amplitudes 
and occupation numbers for boson states of the form ( 38 1 . 



The pairing approximation for bosons is further ce- 
mented by the observation that, fully analogously to the 



fermion theory, in our mean-field theory for bosons the 
left-hand side of (42) is a constant of the motion. This 



also gives an interesting piece of insight into the prob- 
lem of finding the steady state: There is a large (infinite) 
number of constants of the motion, so that time evolution 
cannot lead to a unique steady state without an explicit 
specification of the values of the constants. Physically, on 
the other hand, it is generally up to the interactions with 
the environment to force specific values, such as the zero 
on the right-hand side of Eq. ( 42 1 . In the corresponding 
fermion problem the zero would follow even without the 
interactions with the environment if the system started 
out as a condensate of molecules, and the same holds for 
the boson system. As advertised, our stationary solution 
corresponds to the assumption that all noncondensate 
atoms are correlated as if they came from dissociation of 
molecules. 

Given that A{e) is real and that it must be positive in 
the hmit e ^ oo by virtue of Eq. (l35|, Eqs. (ISSj) and (|42l 



may be solved for the pairing amplitudes and occupation 
numbers, 



Aie) = 
P{e) = 



I3fl 



2v/2(e-/i)2-/32f]2 ' 

1 / / /32f]2 

2 [y 2{e~fi)^-l3^n^ 



1 - 1 



(43) 
(44) 



Hence, given the chemical potential and the ampli- 
tude of the molecular condensate /3, both the occupation 
numbers P(e) and pairing amplitudes A(e) are uniquely 
determined. The inequality 



(45) 
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must hold, else complex occupation numbers would re- 
sult. The equality in Eq. (45 1 presents no problem, since 



the ensuing singularities in the occupation numbers and 
pairing amplitudes are sufficiently mild not to hamper 
the analysis. 



E. Renormalization 



to replace the detuning 6 on the left-hand side with the 
renormalized detuning S. Here we have played fast and 
loose with mathematical rigor, but this could be reme- 
died by introducing the upper limit M to the integration 
just as in Eq. ( 15 1 and then letting M ^ oo. 



Given the occupation numbers ( 44 1 , the integral in the 



normalization equation (30 1 converges as written, but not 



Equations ([33]), (34 1 and (30 1 suffice to determine the Eq. (46) replaces P(e) with 



so the integral involving P{e) in the expression for the 
energy per particle (31). However, if one analogously to 



remaining unknowns a, (3 and ii, although a few issues 
remain. Next we discuss renormalization. 



Consideration of the form of A{e) in Eq. (43 1 shows 
right away that the integral / dey/eA{e) in Eq. (34 1 di- 
verges. However, it turns out that the same renormaliza- 
tion that we devised for the two-atom case also resolves 
this divergence. We replace A{e) with 



Pie) = P{e) 



my 

8e2 



(47) 



in the integral involving A{e) does the subtraction (46 1, 
and replaces the detuning with the renormalized value S, 
all divergences in Eq. (31 1 cancel. We regard this as an 



A{e) = 



(3n 



139. 
2'J2t 



(46) 



which makes the integral convergent. But, to keep 
Eq. (34 1 valid, we need to add the the divergent integral 
/p°°dee~^/^ with an appropriate factor to the left-hand 
side as well. It turns out that the net effect is precisely 



impressive demonstration of the consistency of the mean- 
field theory. 



F. Statement of numerical problem 



We now have deal with the equations 



< a ' 



a H TTTT^Al/ 



Aij2(jn) 
Pi/2im) : 
P3/2{m) -■ 



V2 [ 423/4e 
dx 




-m — 2x171 + 1 



dx 



dx 



^/x^/[{m + x)'^ - 1] [x + ^{m + xf - 1 
\[x_ , 

^ (to + xY' — 1 ^TO + .T + ^ (to + xY' — 1^ 

— Stox^ — 4TO^a::^ + 3a:^ — 2TOa; — rr? + \ 



2~JxyJ{m^ xY - 1 2(a; + TO)a;2 + (2x2 + 1)^(^^7+2^)2^^^ 



^1; (48) 
(49) 
(50) 

(51) 
(52) 



Equations (48 1 are the ones to solve for the unknown 



III. SOLVING THE THEORY 



quantities a, /? and and Eq. (49 1 gives the resulting 
mean-field energy per particle. The integral Aij^ijn) is 
a representation of the integral / dt e^/2^(g) as a func- 
tion of the chemical potential /i, and similarly for P1/2, 
P3/2- These integrals are properly renormalized and di- 
mensionless, and we have gone so far as to write them 
in forms that do not involve near-canceling subtractions 
of large numbers. They are suitable for use in numerical 
computations as written. 



We have solved Eqs. (48l-(52l using Mathematica 



in a combination of analytical and numerical calcula- 
tions. We double-checked many of the results by indepen- 
dent programming on Maple |42j . Unlike in the fermion 
case |31j where the production of accurate numerical re- 
sults for arbitrary parameter values was a major project 
in numerical analysis, with bosons we never had to re- 
sort to a general-purpose programming language such as 
C++. 



8 



In detailed studies we first investigated numerically 
how the results behave, and use this knowledge to for- 
mulate Ansatz solutions to find analytical results. The 
discussion preceding Eqs. (60 1 below serves as an exam- 
ple. An enormous amount of detail could be extracted 
in this way, but our aim is to demonstrate a few major 
qualitative features only. 



Atomic condensate present 



One way of satisfying Eq. ( 33 1 , or the same equation 
as the first member of Eqs. (48), is to require that the 



(53) 



expression inside the brackets vanishes. This leads to 

the integrals become ^1/2(1) -2%/2, Pi/2(1) = 2^2/2, 
and -P3/2(l) = — 8\/2/5, and the equations to solve for a 
and P are 



a 



3/33/2^(5/2 



27/4 



3/2 



3/2 



25/44/2 



1. 



(54) 



(55) 



In principle these have an explicit closed-form solution, 
but in practice we have found it useless and proceed nu- 
merically. Given the solution, the energy per particle 
may be found from 



3(-/i) 



5/2 



5V2 



,3/2 



(56) 



Now, by eliminating a from Eqs. (54 1 and (55 1 we have 



a necessary condition for the solution. 



0- 



,2 a/2 {V2f3n - 



S)/3 23/4(/3r!)3/2 



3/2 



(57) 



It is easy to see numerically that for any real S and > 
this has at most one solution with < /3 < 1. Moreover, 
numerically one may demonstrate that such a solution 
/9(f2, S) is a decreasing function of 6. Therefore, a lower 
limit on the existence of a solution as a function of S is 
a possibility: Once (3 has reached a certain value < 1, 
Eq. (55 1 dictates that = 0. Suppose one decreases 5 



further, then /3 would increase further and Eq. (55 1 would 



require < 0, which is not allowed for a real a. 

The limiting case is found setting a = in Eqs. (53 1- 
( [55] ), which gives a relation between S and fi. We plot in 
Fig. [2] the detuning 6 obtained numerically in this way 
as a function of the Rabi frequency fi, effectively using 
the Fermi frequency as a scale of frequencies. Given the 
equations to solve, we may also attempt power series so- 
lutions. For instance, we may insert the series 

-5 = Kifl^ + K3/2fl^/'^ + K2fl'^ 



Af3/2f^3/2 



(58) 
(59) 



to 




FIG. 2: The line in the {Q, 5) plane separating the phases of 
the system with a BEG of atoms present (a 7^ 0) and absent 
(a = 0). 



into Eqs. (53l-(55), and attempt to satisfy the equations 



power by power in f2. This Ansatz works and gives solv- 
able equations for the coefficients K and L. Following 
such principles, we find the expressions of the curve in 
Fig. |2] for the limits of both small and large fi. 



2. Atomic condensate absent 



(60) 



The other way to satisfy Eq. ( 33 1 is to set 



a = 0, 



(61) 



which corresponds to the up-front statement that there 
is no atomic condensate. Equations (48) then turn into 



{S-2^i)/3 



213/44/2 ' 

3(/3fl)3/2Pi(TO) 



/3 



211/44/2 



(62) 
(63) 

(64) 



which are to be solved for (3 and /i; we have expressed 
the chemical potential /i using the dimensionless variable 
m £ [1,00). The integrals (50)-(52l may be written in 



terms of elliptic integrals, but this fact appears to be 
useless and the practical solutions again proceed numeri- 
cally. Once the solution is found, the energy per particle 
is given by Eq. (49) as written. 



Now, consider Eqs. (62|-(64l for a fixed value of the 



parameter f2, regarding /3 and S as functions of the vari- 
able m. By plotting the respective functions, it may 
be seen that, for m > 1, Pi/2(to) > 0, P{, 2(771) < 0, 
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^i/2('Ti) < 0, and A'^i^{m) < 0, with ' denoting the 
derivative. Equation (64 1 then imphes that (3'{m) > 0, 
and Eq. (63 1 consequently gives S'{m) < 0. In other 



words, when the parameter S is increased while keeping 
fixed, m and f3 resulting from Eqs. (62)~(64l decrease. 



But, by the time the detuning S has reached the line in 
Fig. [2] the parameter m has attained the minimum per- 
missible value TO = 1, and there cannot be a solution for 
any larger S. 



3. Role of atomic condensate 

The two classes of solutions we have found, a > and 
a = corresponding to the presence and absence of an 
atomic condensate join continuously. On the line drawn 
in Fig. [2] a = and /_* = —f3Cl/\^ both hold true, so 
that all of the Eqs. (|53l)-(l55l) as well as Eqs. (|62l)-(l64l) 



are satisfied simultaneously. 

Summarizing, we have the following observations 
about the structure of the theory. There are two dif- 
ferent kind of solutions for fi, a and (3 when the param- 
eters n and 6 are varied, characterized by the conditions 
H = — /3ri/-\/2 and a = 0. Only the former (latter) ex- 
ists in the region of parameters and 6 labeled a ^ 
{a = 0) in Fig. |2] On the borderline between the regions 
both solutions exist and agree, so that they go continu- 
ously from one to the other as the variable 6 and/or 
crosses the line. The solution is unique for all S and 
under the assumptions O > 0, a > and /3 > 0, which 
are a matter of convenience and can always be made. Fi- 
nally, our numerical computations demonstrate that the 
unique solution always exists. 



IV. FEATURES OF THE THEORY 

The present boson theory has three parameters with 
the dimension of frequency, fi, S, and ep- In our dis- 
cussions we regard Fermi energy, , as the scale of fre- 
quencies, although we always write it down explicitly. In 
present-day dilute quantum degenerate gases the repre- 
sentative value is ep ^ ^tt x 10 kHz. In the classic Fesh- 
bach resonance experiments with bosons orders of mag- 
nitude of the coupling Rabi frequencies f2 ~ lOe^ are 
typical |37j . but a narrower Feshbach resonance could 
change this comparison significantly. In photoassocia- 
tion the ratio fl/ep may be varied by varying the inten- 
sity (intensities) of the laser(s). In a Feshbach resonance 
the detuning S depends on the product of magnetic field 
and the difference of the magnetic dipole moments of 
the bound molecular state and of the free-atom state. 
Roughly, 5 ^ ep corresponds to 10 mG change in the 
magnetic field. In photoassociation the parameter 6 is 
directly a matter of tuning of the laser(s). 

For comparison it is useful to recall a few basic facts 
about the usual one-channel theory for a BEG [Ml 112 • In 
such modeling there are no explicit molecules at all, but 




5 S/£p 



FIG. 3: (Color online). Fraction of atoms in the atomic BEC, 
Qf^, as a function of the coupling strength, Q,, and detuning, 
S. 



it is assumed that the atoms interact among themselves 
as characterized by a scattering length a that may be 
tuned by varying the magnetic field around the Feshbach 
resonance, i.e., by varying the detuning. At positive de- 
tunings the scattering length is negative, which indicates 
an attractive interaction between the atoms and collapse 
of the (untrapped) condensate. At negative detunings, 
at least close to the resonance, the scattering length is 
positive and large. According to the theory of molecu- 
lar structure, this means that there must be a weakly 
bound molecular state, which, however, is not included 
explicitly in the one-channel theory. By elementary ther- 
modynamics, at low temperatures these molecules should 
make the thermal-equilibrium state. Positive (negative) 
detunings thus indicate the atom (molecule) side of the 
resonance. 



A. Scaling with detuning and coupling strength 

In Fig. |3] we plot the fraction of atoms in the atomic 
BEG, a^, as a function of the coupling strength and 
detuning S. At f2 <C ef, the system switches abruptly 
(compared to e^) from having all atoms to no atoms in 
the condensate when the detuning crosses the two-body 
resonance position 5 = 0. As the coupling f2 increases, 
the transition from ~ 1 to ^ 1 rounds out. How- 
ever, = is always reached on the line shown in Fig. [2] 
Given that on this line the solution of the system is not 
analytic, we predict a phase transition. 

The analogous plot for the chemical potential /i as a 
function of coupling strength and detuning is given in 
Fig. [4] The function fj,{D,,S) is also nonanalytic on the 
line shown in Fig. |2] but the kink is not visible in Fig. |4] 
The rounding-out with increasing coupling strength is 
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Q/£p 3 



FIG. 4: (Color online). Chemical potential ^ as a function of 
the coupling strength and detuning 5. 



3.0 
2.5 
2.0 
1.5 
1.0 
0.5 
0.0 



I 

1 






\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 

\ 


>p 
<p 

> a 
< a 




~ \ 






■V\ 1 


r ^ — t 1 — 


— 1 — 



FIG. 5: Scaling functions f^{x), etc., in Eqs. (65 1 and (66 1. 
The legend identifies the combination of < or > and a or P 
for each curve. 



again obvious. 

We quantify the rounding-out starting with the numer- 
ical observation that in the regime a ^ Q, the following 
limits hold true: 



J a 

/> 



6 

u 



f. 



(65) 



(66) 



Knowing these scalings it is easy to work out analytical 
expressions for the functions f^, etc., from Eqs. (62 1- 
(64 1; simply put in the scalings, and see what becomes 



of the equations in the corresponding limit of f2. For 
instance, we have fp{x) = [j/(a;)]^, where y is & solution 



to the polynomial equation 

Sy^ - 4^2 x^y^ + 8xy - 2^2 ^ (67) 

for the fixed value of x. We have plotted the scaling 
functions in Fig.jsj They tend to the proper limits in both 
ends of the range of the variable x, and, as appropriate for 
useful dimensionless scaling functions, their magnitudes 
and scales of variation with the argument x are both on 
the order of unity. 

We next turn to the behavior of the chemical potential 
and the energy per particle, particularly in the region a — 
in Fig.|2]with both ^ < and ^ < 0. To this end we first 



note the asymptotic expansions of the functions ( 50 1-( 52 1 
for 771 — > oo: 



^l/2(™) 

Pi/2{m) 
P3/2{m) 



^ 16 



-3/2 + 0(to-7/2), (68) 
5/2 + 0(^-9/2), (69) 



.^^+^„,-V^ + 0{m-y^). (70) 
4 ^ 128 ^ ' y ' 



4x/r7i 128 



The coefficients were inferred from numerical results as 
given by Mathematica's Nlntegrate function. They are 
probably exact; the coefficients of the leading terms ap- 
pear to be precise to at least 10"^, and the coefficients 
in the next-to-leading terms to about 10^"^. 

Now, inserting power series expansions of the form 



2 



^ = _ + Mi/2|^|i/2^A/o|^r- 



/3 = 1 + S_i/2|^| 



-1/2 



B^i\~S\-' 



(71) 
(72) 



into Eqs. (62)-(64l and e xpan din g the result to a power 
series in lj\d\ using Eqs. (68h(70|, we discover that we 
may determine the coefficients M and B in Eqs. (71 1 

_3 and -B_4 in terms of 
([68|-(|70l). These, in turn. 



and ( 72 ) down to and including 
icit coefficients in Eqs 



the exp 

will give an expression for the energy per particle ( 49 1 in 
the form 



+ i?i/2|^r/2 + £;o|^|° + 



(73) 



also down to and including the order \6\ 
We will not write down the expansions 



71 , (721 



and (73 1 in detail, but a few notes are relevant. The ex- 



pansions of chemical potential and energy go with powers 
of the dimensionless number z — n'^/ep\S\ down to and 
including the order and the expansion of conden- 

sate amplitude down to and including the order |i5|~2, 
although in the expansions of fi and e the coefficients of 
\S\~^ happen to be zero, M_i = E^i = 0. But since we 
have fl oc pi/2 = (7V/y)i/2 and ep oc p2/3 ^ (N/V^/^, 
the parameter z does not depend on density, or particle 
number, or volume, at all. It is a single-molecule quan- 
tity. In fact, with the identification (A9), down to the 
order \S\~^ the expansions of /i and e coincide with the 
expansion of Wb/2, where LOf, is the bound-state energy of 
a single molecule as given in Eq. (A15) in the Appendix. 
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Many-body phenomena start at the order |^|~'^/^, from 
which onwards other combinations of f2 and ep than z 
also emerge. For instance, asymptotically, \S\ — + oo, the 
dependence of chemical potential and energy on detuning 
and density goes as (x p/\6\'^^^. The analogous mean-field 
theory for fermions [31] behaves in this respect in the 
same way. 

The conventional scattering length for atoms behaves 
as a cx — (5~^ near the Feshbach resonance, and for 
fermions the scattering length of the molecules near the 
resonance should be 0.6 times the scattering length for 
the atoms [IT]. This means that the density dependence 
of energy should emerge in the order S~^. The present 
mean-field theory does not conform with the standard 
expectations. Whether this is a contradiction or not, we 
cannot say. First, by their very nature the asymptotic 
expansions are valid far away from the resonance, not 
close where one expects the 0.6 a for molecules made of 
fermions. Second, as a practical matter, the frequency 
^^/ep is typically quite large for a Feshbach resonance; 
we estimate ^ 27rxl0 GHz for the usual Feshbach reso- 
nance in ^^Rb [15], and indeed THz scale values for the 
834 G Feshbach resonance in fermionic ^Li [31]. Our 
asymptotic expansions only become useful at extremely 
large detunings, when in practice some other physics as- 
sumption of our theory such as the contact interaction 
or the neglect of the background scattering length has 
already become invalid. We do not expect the expan- 
sions (71|-(73l to be of much practical value, which is 
one reason why we have not listed the coefficients. We 
will return to their theoretically interesting properties 
below, though. 

We next turn to the limit of strong coupling, — > oo. 
Let us first 6x S — 0, so that the ar gum e nt is within the 
regime a ^ 0. We set (5 = in Eqs. (53l-(55), attempt a 
solution of the form 

a = Ao + A^2n-^ + A_4f^"* + . . . , (74) 
-/i = Mo + M^2^-^ + M^i^-^ + . . . , (75) 

and find that it works; asymptotically, with — > oo, 
we have ^ = —^ep and = 3/4. From Eq. (56), the 
corresponding expression for energy per particle is e = 
— |eF. As expected on the basis of unitarity [33] H5] . 
in the limit of very strong interactions the interaction 
strength vanishes from the result and the only energy 
scale that remains, ep, is set by the density of the gas. 

Given that the line separating the a = and a =/= 
regions in Fig. [2] may, perhaps, be thought of as the 
position of the Feshbach resonance shifted by many-body 
effects, it is also instructive to find the asymptotic limits 
of chemical potential and energy as one moves along this 
line. To find them, we set a = in Eqs. (53l-(55) and 
put in the Ansatz 



-6 



+ Don" 4 



(76) 
(77) 



Aside from reproducing the second of Eqs. (60), this 



strength we have /i 
e= -3v^/5eF. 



\/2ep. This, in turn, gives 



B. Thermodynamics 

The atom-molecule BEG has some interesting thermo- 
dynamical properties, in particular a propensity for neg- 
ative pressure. We now discuss these aspects. 

It seems natural to identify the mean-field value of the 
energy E — KNe as the thermodynamic internal energy 
U of the system. Moreover, the system is at zero tem- 
perature, which indicates an entropy identically equal to 
zero. In the sense of thermodynamics, the chemical po- 
tential should then equal 



dN L 



(78) 



There is a subtlety to this definition, in that N is the 
invariant atom number that combines both atoms and 
molecules, yet we have used it to find the chemical po- 
tential for atoms only. On the other hand, the condition 
for chemical equilibrium for atoms and molecules dictates 
that the thermodynamic chemical potentials for atoms 
and molecules satisfy 2/iT.a = ^J■T,m■ From this observa- 



tion it is easy to deduce that the procedure (78) gives 
the correct chemical potential for the atoms. There are 
related issues elsewhere in our thermodynamics discus- 
sions that could be resolved similarly, but henceforth we 
will not bring them up. 

We have talked of the quantity ^ as the chemical poten- 
tial, which is consistent with our thermodynamics iden- 
tification E ^ U if 



f d[Ne] \ 



(79) 



As already noted repeatedly, both atom number and 
volume are ingredients in our theory because the param- 
eters and tp depend on the density; for instance. 



dN 



2N' 



(dep\ 
\dNjy 



2ep 
3N 



(80) 



However, verifying Eq. (79) by simply taking the an- 



alytical N derivative of Eq. (49) appears cumbersome. 



analysis shows that in the limit of large interaction 



since we should then take into account the implicit de- 
pendence on N of the quantities a and /3, and even of /i 
itself. Instead, we have carried out the verification nu- 
merically, by calculating the derivative of e numerically 
using Eqs. (80 1, and found that Eq. (79) is indeed satis- 
fied. 

We have also found that the analytical asymptotic ex- 
pansions ( 73 1 and ([tTJ satisfy Eq. ( 79 1 to all orders, down 
to and including for which the expansion coeffi- 

cients are explicitly known. This is a rather impressive 
confirmation of our originally numerical identification of 
the expansion coefficients in Eqs. (68)-(70l. 
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Finally, there is the question of pressure. The ther- 
modynamical expression for pressure in the present case 
with identically zero entropy is 



dU_ 
dV 



-hN 



N 



de 
dV 



(81) 



N 



But, given the identifications we have already made, the 
Gibbs-Duhem relation of thermodynamics reads 



or 



hNe = -pV + hNfj, , 



hp 



(82) 



(83) 



We have calculated the pressure both ways numerically. 



too, and found that Eqs. (81 1 and (83) are consistent 



We have already noted that in the case S — and 
Q — + oo, the limits /i ^\^f and e — > —\(-f hold 
true, which gives p/{hp) —^cp- The pressure becomes 
negative. In a way, this is as expected. Assuming that 
the transition line between the a = and a ^ lines 
denotes the true Feshbach resonance as shifted by the 
interactions in the system, this limiting case is on the 
atom side the Feshbach resonance where, according to 
the standard one-channel model, the scattering length is 
negative and the BEG of the atoms is liable to collapse. 
This is also what one would expect for negative pressure. 

The situation becomes more peculiar after a look at 
Fig.|6] which plots the negative of the pressure for a range 
of the parameters and S. Although the absolute value 
of the pressure drops sharply toward of the borderline 
between the areas of a = and a ^ when entering 
from the a ^ side, as far as we can tell, the pressure 
remains negative for all values of f2 and S. For instance, 
the expansions ( 73 1 and ( 71 1 give 



P_ 

hp 



0{\6\-'). 



(84) 



For comparison, we have plotted in FigJTlthe pressure 
of a two-component Fermi gas from Eq. ( |83[ ) in complete 
analogy with Fig. |6] No negative pressure develops for 
a Fermi gas in the corresponding mean-field theory |31j . 
Moreover, there is no question about the experimental 
stability of the Fermi gas in the vicinity of the Feshbach 
resonance. 

As it comes to the Bose gas, the negative pressure is 
a mixed blessing. It is an interesting prediction in its 
own right, but complicates the observation of the other 
features we have discussed. We hope that additional 
atom-atom, atom- molecule, and molecule-molecule inter- 
actions not included in our model might stabilize the 
Bose gas sufficiently, especially on the a — side, that 
the phase transition like feature and the many-body shift 
of the Feshbach resonance could be seen under favorable 
circumstances. 
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FIG. 6: (Color online). Pressure p of the atom-molecule BEC 
plotted as a function of the interaction strength fl and detun- 
ing 5. 
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FIG. 7: (Color online). Pressure p of two-component Fermi 
gas plotted as a function of the interaction strength Q and 
detuning 5. The definition of the interaction strength SI is 
as in Ref. [31], but, for compatibility with the present paper, 
the detuning here is the detuning of Ref. .31 divided by two. 
For comparison, for an ideal Fermi gas at zero temperature, 
P ■- 



V. CONCLUDING REMARKS 

So far the molecular lifetimes have not allowed atom- 
molecule equilibration in a Bose gas near a Feshbach 
resonance, so that the kind of systems we have stud- 
ied are presently not feasible experimentally. However, 
progress is continually made j47j in two-photon (two- 
color) photoassociation in a Raman scheme [B], and in 
heteronuclear systems direct one-photon photoassocia- 
tion |48j from the dissociation continuum to a low-lying 
vibrational level is possible in principle. We are of the 
opinion that photoassociation to a stable molecular state 
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will eventually be achieved in a BEC. 

There is a semantic issue with the mean-field theory, in 
that, for instance, the review [2U] regards our approach 
as going beyond the mean-field theory. Now, in a usual 
atomic BEC "mean-field theory" means Gross-Pitaevskii 
equation for the condensate [33]. The analogy in atom- 
molecule systems would, indeed, be an analysis in terms 
of atomic and molecular condensates only, with no non- 
condensate atoms or molecules. This was the case in the 
early discussions jl7lll9j . which we nowadays term "two- 
mode" models. 

However, such an analysis ignores an inherent asymme- 
try in atom-molecule conversion. According to momen- 
tum conservation, in a Feshbach resonance in free space 
two zero-momentum condensate atoms may only combine 
into a zero-momentum condensate molecule, but in the 
reverse process the two atoms dissociated from a zero- 
momentum molecule may have arbitrary opposite mo- 
menta. Taking into account such "rogue" dissociation [37| 
was the original aim of our present modeling that ac- 
counts for correlated pairs of atoms. Initially we regarded 
this as an extension of mean-field theory, but after we dis- 
covered that the time independent version for fermions is 
a variation of the BCS theory, one of the quintessential 
mean-field theories, we refer to our approach simply as 
mean-field theory. 

Mean-field theories, of course, are an interesting case 
in their own right. They are typically the first, and oc- 
casionally the last, tool in the analysis of a new phase 
transition, even though it is well known that they can- 
not be expected to be quantitatively accurate. In fact, 
we cannot think of any instance in the area of quantum 
degenerate dilute gases in which a properly formulated 
mean-field theory has given a qualitatively wrong pre- 
diction, and only very recently has it become possible 
to distinguish between mean-field theory and a strongly 
correlated approach quantitatively in an experiment; c.f. 
Refs. [H] and [IS]. 

These observations cast an interesting light on our pre- 
diction of negative pressure even on the molecule side 
of the Feshbach resonance, where the standard picture 
is that the scattering length should be large and pos- 
itive. Granted, our calculations are missing the usual 
background scattering length and are cursory about the 
structure and collisions of the molecules. One may also 
argue that, as we have no time scale to offer for the as- 
sociated instability, it may be unobservable even if it re- 
ally existed. After all, at sub-Kelvin temperatures, at 
any pressure, the thermodynamic ground state of alkali 
metals is a solid not a gas, which has not precluded in- 
numerable successful BEC experiments. 

Nonetheless, the negative pressure may be viewed as 
an opportunity rather than a nuisance. It either exists, 
or is a qualitatively wrong prediction from a mean-field 
theory. Speculating further, since we do not have any 
time scale for the collapse associated with the negative 
pressure, we do not know that it is short either. It is 
a well-known empirical fact that in the neighborhood of 



a Feshbach resonance, on both sides of the resonance, 
atoms are lost from a BEC. Perhaps the negative pressure 
contributes to, or even dominates, the loss. 

We close with a remark about a rather esoteric original 
motivation for this particular piece of work on bosons. 
Some time ago we noticed |49l in a simple model for 
two trapped ions that the statistics has no effect on the 
thermodynamic properties. Our interpretation was that 
the Coulomb interaction between the ions keeps them 
sufficiently far apart, so that they are effectively distin- 
guishable and quantum statistics is moot. Calculations 
in which strongly interacting bosons crystallize |50j in a 
trap just like ions do then led us to the idea of "superuni- 
versality:" If strong interactions between the particles 
keep them apart, the state could not only not depend on 
the strength of the interactions, but could also be inde- 
pendent of atom statistics. Now, within the confines of 
mean-field theory, our results for bosons and fermions are 
quite different. Technically, superuniversality in mean- 
field theories is neither a necessary nor a sufficient condi- 
tion for superuniversality in Nature, but nonetheless, the 
notion of superuniversality did not pass this particular 
test. 
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APPENDIX A: STEADY STATE 

In this Appendix we study the steady state of the two- 
atom model as a bound-continuum problem along the 
lines of the classic paper by Fano [40 , with due consid- 
eration to renormalization of the ultraviolet divergence. 

We imagine that there is initially a bound molecu- 
lar state 1 6) and a continuum of states |e) labeled by 
frequency and normalized in such a way that (e|e') = 
6{e— e'), coupled by a photoassociative or Feshbach res- 
onance coupling. We write the Hamiltonian as 



H 
with 



\b)S(b\+ / de\e)S(e 



K{e) 



deK{e) [\b){e\ 



k>(6|] , 

(Al) 
(A2) 



The relevant features in Eq. (A2) arc the fourth root of 



frequency (energy), and the constant k with the dimen- 
sion of frequency; the rest of the constants are an attempt 
to simplify the appearance of a few results below. Simi- 
larly, the state vector is 



= / dea(e;i)|e). 



(A3) 
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The time dependent Schrodinger equation becomes 
. d 



+ / deK{e)aie), 



d 

z — a(e) = ea{e) + K{e)h, 



(A4) 
(A5) 



and the time independent Schrodinger equation for an 
eigenstate of the Hamiltonian with the eigenfrequency lj 
is 



(A6) 
(A7) 



{uj~5)h = J deK{e)a{e), 
(w-e)a(e) = K{e)b. 



The substitution 



a(e) 



.3/2 



A 



(i) 



(A8) 



converts Eqs. (A6l and (A7) to Eqs. (12 1 and dl3 



Moreover, let us solve Eqs. (A4| and (A5) in the pole 
approximation for a positive detuning 6, and likewise 
solve Eqs. (10 1 and (111 in the pole (Wigner-Weisskopf) 



approximation, then the ensuing decay rates for the 
molecules are equal if the parameters are related by 



5124 



(A9) 



The present continuum problem ( Al I simply solves the 



dilute-gas (single molecule) limit of our mean-field the- 
ory, but with the advantage that the amplitudes b and 
a(e) also have the normalization conditions that follow 
from the orthonormality conditions of the states \b) and 
|e) and the manifest hermiticity of Eqs. (A4| and ( |A5[ ). 
In the steady state the amplitudes b and a(e) , of course, 
depend on the eigenvalue lu, a dependence that we will 
write down below. 



Now, the formal solution for a(e) of Eq. (A7l, 



a(e, ll>) 



K{e)b{Lo) 



(AlO) 



with 9 being the usual unit step function. The integral on 
the right-hand side has the same ultraviolet divergence 
as before, and the cure is exactly the same; we add the 
infinity / deK^{e)/e to both sides of the equation, which 
renormalizes the detuning and makes the principal value 
integral convergent. 



to - S = d{Lo)f{uj)K^{uj) +V 



e(w - e) 



1. Bound State 

Suppose first that w < 0. Then the principal value 
integral is a usual integral, and we have the equation 



(A14) 



It turns out that this equation has a real solution if and 
only if (5 < 0, and the unique solution is then 



UJb 



= S -2k + 2\J - k6 . 



(A15) 



There is one, and only one, negative-energy solution if 
and only if the detuning is negative. 

One might surmise that the negative-energy solution 
is bounded, i.e., normalizable to unity. In fact, using 
Eq. (TAToI), we find 



de \a{e,ujb)\" 



-UJb 



(AI6) 



so that the bound and continuum amplitudes in the wave 
function normalized to unity are 




b{ujb) 



a(e,u;b) = b(uib) ■ 

UJb — e 



2. Continuum States 



(A17) 
(A18) 



contains a singularity at e — uj that renders the mean- 



ing of the right-hand side of Eq. (A6I ambiguous. The 
substitution u uj — iri, with 77 = 0-I-, removes the ambi- 
guity, but gives at most one stationary state. Following 
Fano j40j, we therefore attempt a solution in the form 



a(e,w) 



V- 



1 



K{e)b{Lo), (All) 



where V denotes the principal value integral and /(w) is 
yet to be determined; f{uj) = in would give the forward- 
in-time solutions we have discussed earlier. Inserting 



Eq. (All I into Eq. (A6) gives 



uj-6^ 0{uj)f{uj)K^{uj) + V de 



(A12) 



The system also has a bountiful of positive-energy 
eigenstates. It turns out that for uj > 



V I de , ^ \ = , 
e{uj - e) 



so that Eq. (A13) gives 



UJ — 6 



(A19) 



(A20) 



For positive energies we are obviously dealing with con- 
tinuum states, so that we aim at the normalization 

b* {uj)b{uj') + [ dea* (e, uj)a{e, uj') = S{uj - uj') . (A21) 
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Note that while for a unit-normalized bound state the 
coefficient b{ujg) is dimensionless, here the dimension of 
b{uj) will be the inverse of the square root of frequency. 
The calculation based on Eq. ( All ) is straightforward ex- 



cept for the ensuing product of principal- value integrals, 
which may be handled as in |40j : 



1 



1 



1 



CO — e uj' — e 

2 



1 



-V- 



1 



+'n:^6{uj - J)5{e - uj) 



(A22) 



The principal value integrals arising in this way give zero 
just as in Eq. (A19l, so we have 



1 . 



(A23) 



The continuum state vectors for lo > are therefore fully specified by 



— (5) + 4^0; 



6M, 



(A24) 
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